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ABSTRACT 

We used the Australia Telescope Compact Array to map a large field of approximately 
2° x 2° around the Sculptor group galaxy NGC 300 in the 21-cm line emission of neutral 
hydrogen. We achieved a 5a H I column density sensitivity of 10 19 cm -2 over a spectral 
channel width of 8 kms -1 for emission filling the 180" x 88" synthesised beam. The 
corresponding Hi mass sensitivity is 1.2 x 10 5 M Q , assuming a distance of 1.9 Mpc. 
For the first time, the vast H I disc of NGC 300 has been mapped over its entire extent 
at a moderately high spatial resolution of about 1 kpc. 

NGC 300 is characterised by a dense inner H I disc, well aligned with the optical 
disc of 290° orientation angle, and an extended outer H I disc with a major axis of 
more than 1° on the sky (equivalent to a diameter of about 35 kpc) and a different 
orientation angle of 332°. A significant fraction (about 43 per cent) of the total detected 
H I mass of 1.5 x 10 9 M Q resides within the extended outer disc. We fitted a tilted ring 
model to the velocity field of NGC 300 to derive the rotation curve out to a radius of 
18.4 kpc, almost twice the range of previous rotation curve studies. The rotation curve 
rises to a maximum velocity of almost 100 kms" 1 and then gently decreases again in 
the outer disc beyond a radius of about 10 kpc. Mass models fitted to the derived 
rotation curve yield good fits for Burkert and NFW dark matter halo models, whereas 
pseudo-isothermal halo models and MOND-based models both struggle to cope with 
the declining rotation curve. 

We also observe significant asymmetries in the outer Hi disc of NGC 300, in 
particular near the edge of the disc, which are possibly due to ram pressure stripping 
of gas by the intergalactic medium (IGM) of the Sculptor group. Our estimates show 
that ram pressure stripping can occur under reasonable assumptions on the density 
of the IGM and the relative velocity of NGC 300. The asymmetries in the gas disc 
suggest a proper motion of NGC 300 toward the south-east. At the same time, our 
data exclude IGM densities of significantly higher than 10~ 5 cm' 3 in the vicinity of 
NGC 300, as otherwise the outer gas disc would have been stripped. 

Key words: galaxies: individual: NGC 300 - galaxies: kinematics and dynamics - 
galaxies: structure - radio lines: galaxies. 



1 INTRODUCTION 

At distances in the range of about 2 to 5 Mpc, the Sculptor 
group i s among the nearest galaxy groups beyond the Local 
Group jjerien. Freeman fc Bi ngacli 1993). It forms an elon- 
gated filament of galaxies and comprises a number of sepa- 
rate subgroups at different distances along the line of sight. 
At a distance of about 2 Mpc, the nearest of these subgroups 
consists of NGC 55, NGC 300, and p ossibly two or more 
know n dwarf spheroidal companions (|Karachentsev et al.l 
120031 : Koribalski et al., in prep.). NGC 55 and NGC 300 are 
medium-size d spiral galaxies of type SB (s)m and SA(s)d, 
respectively (|de Vaucouleurs et al.l Il99lf ) . Their proximity 
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makes them preferential targets for deep Hi observations 
with high spatial resolution and Hi mass sensitivity. 

Some of the basic properties and physical parame- 
ters of NGC 300, including the results of this work, are 
listed in Table Q] The distance towards NGC 300 has 
been mea sured with g reat a ccuracy and through different 
methods. iRizzi et alj ()2006h used the tip of the red gi- 
ant branch to determine distance moduli of (m — M)o = 
26.30 ± 0.03 ± 0.12 mag and 26.36 ± 0.02 ± 0.12 mag 
through two differ ent statistical metho ds as part of the 
Araucaria Project (|Gieren et alj l2005bl ). Their results are 
consistent with the distance moduli of 26.43±0.04±0.05 and 
26.37 ± 0.05 ± 0.03 mag derived from Cephei d variables by 
iGieren et all l|2004l) and iGieren et all i|2005al h respectively. 
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Table 1. Properties of NGC 300. The separation between in- 
ner and outer Hi disc was made at a column density level 
of 5 X 10 20 cm -2 . Note that all errors are statistical uncer- 
tainties and do not reflect possible syst ematic errors . Refer - 
en ces: [Hide Vaucouleurs et all Jl99lt); [2 1 iGieren et all (|2005eJ) : 
[3] iRizzi et al. I <2006l) ; [41 ICarignanl [I1985T) . Values without refer- 
ence are the result of this study. 



Parameter 


Value 


Unit 


Ref. 


iype 


S A 

O/i^a yo 




Ml 


a (J2000) 


00 h 54 m 53?4 






S (J2000) 


-37°41'02'.'6 






Distance 


1.9 


Mpc 


[2] [3] 


R2S 


9'.8 




[4] 


Radial velocity 








baryccntric 


144 ±2 


kms -1 




LSR 


136 ± 2 


kms -1 




Position angle 1 








inner disc 


290?0±0?5 






outer disc 


331?7±0?3 






Inclination 


40° . . . 50° 






Max. rot. velocity 


98.8 ±3.1 


kms -1 




Integrated flux 


1.72 x 10 3 


Jy km s 




H I mass 


1.5 x 10 9 


M 




inner disc 


8.5 x 10 s 


M 




outer disc 


6.5 x 10 s 


M 




Gas mass 2 


(1.9 ±0.2) x 10 9 


M 




Stellar mass 2 


(1.0 ±0.1) x 10 9 


M 




Total mass 3 


(2.9 ±0.2) x 10 10 


M 





1 w.r.t. the J2000.0 equatorial coordinate system 

2 across our tilted ring area (0.9 < r < 18.4 kpc) 

3 within r = 18.4 kpc 



Based on these measurements we therefore adopt a distance 
of 1.9 Mpc for NGC 300 throughout this paper. 

Early H I imaging of N GC 300 was carried out by 
IShobbrook fc Robinson! <| 19671 ) using the 64-m Parkes radio 
telescope. They determined an Hi mass of 2.1 x 10 9 M 
and a total mass of (2.5 ± 0.5) x 10 10 M from the rota- 
tion curve under the assumption of a distance of 1.9 Mpc. 
They also noticed that the velocity field of NGC 300 appears 
strongly distorted which they attributed to the gravitational 
influence of a large H I cloud to the south-east of the galaxy. 
This cloud, however, w as later found to be most likely part of 
the Magellanic Stream |Mathewson. Clearv fc Murravl 19751 ; 
IRogstad, Chu fc Crutcherlll979l '). ~ 

A larger area aroun d NGC 300 was mapped in Hi by 
iMathewson et all (|l975l ) with the 64-m Parkes radio tele- 
scope. Within their brightness temperature sensitivity of 
0.3 K over 4.1 kms -1 and across the 15 arcmin Parkes beam 
they detected several extended gas clouds with velocities 
comparable to those found across NGC 300, including a long 
'H I ta il' extending abou t 2° to the southeast of the galaxy. 
While IMathewson et all (|l975l ) assumed that these clouds 
were associated with NGC 300, there has been a long debate 
over the past decades whether they might instead be frag- 
ments of the Magell anic Stream that runs across the same 
part of the sky (e.g. lHavnes fc Robertsiri979l ; IPutman et alj 
l2003h . 

A more detailed H I image of NGC 300 was obtained by 
IRogstad et ail (|l979h with the twin-element interferometer 
of the Owens Valley Radio Observatory. At a much higher 
angular resolution of 2' x 3' half-power beam width they were 



able to derive a detailed radial velocity map of the inner part 
of the galaxy. By fitting a tilted ring model to the velocity 
field they were able to describe the distorted velocity con- 
tours of NGC 300 as a result of varying inclination and posi- 
tion angle acro s s the H I disc. The rotation curve obtained by 
IRogstad et all l|l979h from the tilted ring model is basically 
flat beyond a radius of about 10 arcmin with a turnover ve- 
locity of 94 kms -1 at 16 arcmin radius. They suggest that 
the warping of the Hi disc of NGC 300 could have been 
caused by a close encounter with another galaxy or massive 
H I cloud within the past 10 9 years. Similar H I synthesis ob- 
servations of NGC 300 with somewhat higher angular reso- 
lution (50" x 50" half-powe r beam width) were obtained by 
IPuche. Carignan fc Bosnia! fl990) who employed a mosaic 
of five pointings with the Very Large Array (VLA). 

Unfortuna t ely, the synthesis images obtained by 
IRogstad et"all l| 19791 ) and (Puche et all l|l990t ) did not re- 
veal the extent and structure of the outer disc of NGC 300 
due to insufficient field of view and lack of sensitivity. We 
therefore decided to obtain deep H I observations of a large 
field of 2° x 2° around NGC 300 in a mosaic of 32 pointings 
with the Australia Telescope Compact Array (ATCA). The 
aim of these observations was to map the entire extent of 
the Hi disc of NGC 300, determine its structure, and dy- 
namics, and search for extra-planar gas in the vicinity of 
the galaxy. The observations and results of this project are 
described in this paper. Similar observations of the Sculptor 
group galaxy NGC 55 will be presented in a separate paper. 

This paper is organized as follows: Section [2] describes 
our observations and data reduction procedure. In Section [3] 
we discuss our results, including the general physical pa- 
rameters of NGC 300 and the results of the rotation curve 
analysis. In Section [4] we describe the fitting of mass mod- 
els, with and without dark matter, to the observed rotation 
curve. In Section [5] we discuss the possible origin of the ex- 
tended outer gas disc of NGC 300 and present evidence of 
the distortion of the gas disc by ram pressure. Finally, Sec- 
tion [6] gives a summary of our results and conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations of NGC 300 were carried out in late 2007 
and 2008 with the six 22-m antennas of the Australia Tele- 
scope Compact Array. Five of the antennas are moveable 
along an east-west track with a total length of 3 km, whereas 
the sixth antenna is fixed and separated by another 3 km 
from the western end of the track. 

The total observing time of 96 h was equally divided 
among the EW352 and EW367 array configurations to im- 
prove baseline coverage. Because of the large spatial sepa- 
ration between antenna 6 and all other antennas in these 
configurations we only used data from antennas 1 to 5 in 
our analysis. The resulting minimum and maximum base- 
lines were 31 and 367 m, respectively. 

NGC 300 was simultaneo usly observ ed in the 21-cm line 
emission of neutral hydrogen l| Wild! 1952! ) and in 21-cm radio 
continuum emission, using the two independent frequency 
chains of the ATCA. For the Hi observations we used a 
bandwidth of 8 MHz centred at a frequency of 1420 MHz. 
The resulting velocity coverage is about 1600 kms" 1 with 
a channel separation of 3.3 kms -1 and an effective velocity 
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Figure 1. Outline of our ATCA mosaic of 32 pointings around 
NGC 300 (north is up). Pointing centres are marked by black 
points with the dotted circles indicating the half-power beam 
width of 33 arcmin of the ATCA antennas at A = 21 cm. The 
black cross marks the central position of the entire mosaic (and 
of NGC 300) at a = 00 h 54 m 53?5 and 8 = -37°41'04" in J2000 
equatorial coordinates. The grey ellipse outlines the approximate 
size and orientation of the optical disc of NGC 300. 



resolution of about 4 km s _1 . The radio continuum data were 
obtained over a bandwidth of 128 MHz centred at 1384 MHz. 
Only the H I data will be discussed in this paper. 

In order to map a large field of about 2° x 2° (cor- 
responding to a projected size of about 65 x 65 kpc 2 ) we 
covered NGC 300 with a large mosaic of 32 separate point- 
ings. Adjacent pointings were separated by 16' 5 which cor- 
responds to approximately half the half-power beam width 
of the ATCA antennas at a wavelength of A = 21 cm. The 
geometry of the entire mosaic is outlined in Fig. [T] The in- 
tegration time per pointing in each array configuration was 
1.5 h. The resulting spectral rms noise level at 3.3 kms -1 
channel separation is about 5.5 mjy per beam which is equal 
to the theoretical noise level. 

At the beginning of each observing run we spent about 
15 to 20 minutes on the ATCA's standard flux calibration 
source, PKS 1934—638, which was used for bandpass and 
flux calibration of the data. Every 45 minutes during our 
observing runs we integrated for 5 minutes on our gain am- 
plitude and phase calibrator, PKS 0008—421, which has a 
20 cm flux of about 4.5 Jy and a separation from NGC 300 
of about 9? 5. PKS 0008—421 was used to calibrate gain am- 
plitudes and phases as a function of time. 

The data were reduced and analysed w ith the 
Astro nomical Image Processing System (AIPS; iGreisenl 
I1990T I. the Multichannel Image Reconstruction, Im- 
age Analysis and Display ( Miriad) software package 
|Sault. Teuben fe Wrightlll995|'). and the Groningen Image 
Processing System (GIP SY; lAllen. Ekers fe Terlouwl 1 19851 : 
Ivan der Hulst et al.lfl992T ) in a combined approach. We first 
used the task ATLOD to read the raw data into AIPS for 
the purpose of flagging bad pixels in the time-baseline do- 
main using the interactive task SPFLG. Next, we imported 
the flagged data into Miriad. The flagging table produced 



by AIPS was then applied via the Miriad task fgflag 
before proceeding with some additional flagging of shad- 
owed antennas, H I absorption in the primary calibrator, 
PKS 1934—638, and otherwise corrupted visibilities. We 
then carried out the standard data reduction procedure, 
making use of Miriad's dedicated mosaicking capabilities 
which include automatic corrections for primary beam at- 
tenuation. 

For the image cube and beam we applied a robust 
weighting of visibilities with a robustness parameter of zero. 
The resulting synthesised beam is elliptical with a full width 
at half maximum (FWHM) of 180" x 88" and a position an- 
gle of —6° with respect to north. Furthermore, we employed 
a pixel size of 30 arcsec in the image domain and a resolu- 
tion and channel separation of 8 kms" 1 through Hanning 
smoothing in the spectral domain. The resulting spectral 
rms noise level is 3.5 mjy per beam, corresponding to a 5a 
Hi column density sensitivity of 1.0 x 10 19 cm -1 per spec- 
tral channel. This translates into a 5a Hi mass sensitivity 
of 1.2 x 10 5 M at a distance of 1.9 Mpc. 

To reduce the impact of beam sidelobes caused by our 
incomplete aperture synthesis we subsequently deconvolved 
the image data c ube using the CLEAN algorithm origi- 
nally developed bv lHogbom Jl974f). We actually applie d the 
Steer CLEAN algorithm (|Steer. Dewdnev fc Itolll984l ) im- 
plemented in the Miriad task mossdi with a flux density 
cutoff of 10 mjy per beam equivalent to about three times 
the rms noise level at 8 kms -1 channel width. Individual 
channel maps of the final data cube are shown in Fig. [2] 

The zeroth, first, and second moment maps of NGC 300 
were calculated using the AIPS task momnt. For this pur- 
pose, we first multiplied the data cube with the theoretical 
gain across the mosaic to obtain a data cube with a con- 
stant noise level at all positions and frequencies. Next, we 
created a mask from the data cube by Hanning-smoothing 
the frequency axis over three channels, convolving the spa- 
tial dimension with a Gaussian of 5 pixels FWHM, and then 
masking all pixels in the smoothed cube that had a flux den- 
sity value of less than 4.5 mjy per beam. In momnt, only 
those pixels of the original data cube that had not been 
masked in the smoothed cube were used in calculating the 
moment maps. Finally, the moment maps were divided again 
by the gain model to recover the original gain levels across 
the mosaic. 

There are several ways of determining the mean radial 
velocity of an H I line, including the calculation of the first 
moment of the spectrum, the fitting of one or more Gaussian 
functions to the spectral line, or simply by using the velocity 
of the peak of the Hi line. The impact of these and other 
methods on th e outcome of rotatio n curve fits was stud- 
ied in detail bv lde Blok et all (|2008T ) based on a sample of 
galaxies ob served as part of T he Hi Nearby Galaxy Survey 
(THINGS; IWalter et al.ll2008f >. They found the most suit- 
able method to be the fitting o f Gauss-Hermite polynomials 
l|van der Marel fc Franx| [l993) with an additional skewness 
pa rameter, h^, to the sp ectral lines. The fitting function used 
bv lde Blok et all |200cT ) and implemented in the GIPSY task 
xgaufit is a modified Gauss-Hermite polynomial including 
only the zeroth- and third-order terms, namely 



A(y) = ae X p(-y 2 /2) [l + h 3 



2j / 3 -3y s 
\/3 , 



(1) 
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Figure 2. Channel maps of the ATCA H I data cube of NGC 300 in the range of «lsr = 16 to 248 km s~ 1 (corresponding to barycentric 
velocities of 24 to 256 kms -1 ). The greyscale images range from —20 to +100 mjy per beam. The contour levels are —15 (dashed), 15, 
50, 150, and 500 mjy per beam. The increased noise along the edges of the maps is the result of primary beam attenuation. Channels 
containing Galactic H I emission were not deconvolved. 
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Figure 3. Upper-left panel: Optical image of NGC 300 from the Digitized Sky Survey (in logarithmic scaling) with Hi column density 
contours overlaid. Upper-right panel: Hi column density map of NGC 300 derived from the zeroth moment under the assumption 
that the gas is optically thin. The black Hi column density contours in this and the previous map correspond to 0.1, 0.4, 1, 2, 4, 8, 
and 12 X 10 20 cm -2 . Lower-left panel: LSR radial velocity field derived from Gauss-Hermite polynomials. The contours are drawn at 
intervals of 15 kms -1 centred on the systemic velocity of NGC 300 of i>lsr = 136 kms" 1 (bold contour line; equivalent to 144 kms 
in the barycentric reference frame). Lower-right panel: Velocity dispersion map derived from Gauss-Hermite polynomials. The contours 
correspond to 8, 10, 12, 14, 16, 18, and 20 kms" 1 . 



with the substitution y = (x — b)/c. The free parameters of 
the fit are a, b, c, and the aforementioned skewness parame- 
ter, ft.3. The advantage of using Gauss-Hermite polynomials 
instead of pure Gaussian functions is that the former will 
account for asymmetries in the line profile that can be char- 
acterised by /13. In our analysis we used the position of the 
maximum of A(y) to derive the radial velocity field. 



3 RESULTS 

3.1 General parameters 

The total H I column density map of NGC 300, derived from 
the zeroth moment of the data cube, is shown in Fig. [3pl The 

1 Note that the assumption was made that the H I gas is optically 
thin. In the innermost regions of NGC 300 this assumption may 
not be valid, and the true column densities may be higher. 



H I disc appears to be significantly more extended than the 
optical disc and shows a clear separation into an inner and 
outer disc. The inner disc is characterised by high column 
densities of Nm > 5 x 10 20 cm -2 and appears to be aligned 
with the stellar disc. In contrast, the outer disc has lower 
column densities of typically Nm. < 2 X 10 20 cm" 2 and is 
spatially much more extended. Furthermore, the transition 
between the inner and outer disc seems to go along with a 
systematic twist resulting in very different orientation an- 
gles of the two discs. This twist is also evident in the radial 
velocity map in Fig. [3] and the position- velocity diagrams in 
Fig.H 

Owing to our large mapping area and sufficiently high 
sensitivity, this is the first time that the entire extent (down 
to a column density level of about 10 19 cm" 2 ) of the Hi 
disc of NGC 300 has been revealed in such detail. Within the 
10 19 cm" 2 column density level, the major axis of the H I disc 
of NGC 300 spans slightly more than 1° on the sky which 
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Figure 4. Position-velocity diagrams along the major and minor axis of the inner disc (top panels, with position angles of tp = 290° and 
200°, respectively) and along the major and minor axis of the outer disc (bottom panels, with tp = 332° and 242°, respectively). Contour 
levels are 5, 30, 100, 300, and 500 mjy per beam. The dotted line indicates the systemic velocity of NGC 300 of v aya = 136 kms -1 in 
the LSR frame (equivalent to 144 kms -1 in the barycentric reference frame). The offset of corresponds to the centre of NGC 300. 
Emission near i^lsr = kms - is due to Galactic foreground. 



is equivalent to a linear diameter of about 35 kpc. Earlier 
H i ob servations of NGC 300 with the VLA bv lPuche et all 
l| 19901 ) revealed only a small part of the outer disc, mainly 
due to the insufficient size of their mosaic of only five point- 
ings. 

The integrated Hi spectrum of NGC 300 is shown 
in Fig. [5] The total flux derived from our ATCA data 
is 1720 Jykms -1 . The corresponding total Hi mass of 
NGC 300 amounts to Mm = 1.5 x 10 9 M under the as- 
sumption that the gas is optically thin. Since this will not be 
the case in the dense inner regions of the disc, the derived 
mass must be considered a lower limit. In addition, some 
of the diffuse flux of NGC 300 will not have been picked 
up by the ATCA due to the lack of sufficiently short tele- 
scope baselines. The effect of missing short spacings and the 
inability of CLEAN to recover all the flux is visible in the 
channel maps of Fig. [51 with the Hi emission sitting in a 
slightly negative bowl. 

The Hi mass determined from the Hi Parkes All-Sky 
Survey (HIPASS; iBarnes et al.ll200ll) is slightly large r with 
M m = (1.69±0.13) x 10 9 M ffl ((Koribalski et alj|2004 scaled 
to a distance of 1.9 Mpc). This indicates that some flux 
is missing in our ATCA observations, although not much. 
On the other hand, HIPASS suffers from inaccurate band- 
pass calibration resulting in serious artefacts near bright 
and extended sources, and the true Hi mass of NGC 300 
could even be somewha t higher than the one determined by 
Koribalski et all J2004). The tota l Hi mass of 8 x 10 8 Mq 
determined bv lPuche et alJ |l99d ) from their VLA observa- 
tions (again scaled to a distance of 1.9 Mpc) is significantly 



smaller, mainly reflecting the limited size of their mosaic on 
the sky and possibly the lack of short telescope baselines. 

We can also try to assess how the mass is distributed 
between the inner and outer disc of NGC 300. By simply 
defining a column density boundary of iVm = 5 x 10 20 cm -2 
between the two disc components we obtain Hi masses of 
8.5 x 10 8 and 6.5 x 10 8 M Q for the inner and outer disc, 
respectively. Hence, a significant fraction of the Hi mass is 
within the extended outer disc. 

A closer inspection of the column density distribution 
reveals a certain degree of asymmetry in the disc. The south- 
eastern edge of the disc is rather clean and sharp, whereas 
the north-western edge looks more frayed and less well de- 
fined. This is illustrated in Fig. [6] showing the H I column 
density profile across NGC 300 as a function of declination 
offset (i.e. along a line of constant right ascension passing 
through the centre of the galaxy). There is a distinctive 
plateau in the southern part of NGC 300 followed by a sud- 
den and steep drop in column density towards the southern 
edge of the disc. In the northern part, in contrast, we observe 
a steady exponential decline without any obvious plateau or 
drop. 

The velocity dispersion map of NGC 300, derived from 
the fitting of Gauss-Hermite polynomials to the H I spectra, 
is shown in the lower-right panel of Fig. [3] Velocity disper- 
sions across the inner disc are variable and somewhat higher 
than in the outer disc with typical values in the range of 
about 10 to 15 kms -1 . In contrast, the outer disc appears 
homogeneous with dispersions of just under 10 kms -1 in 
most areas (equivalent to just over 20 kms -1 FWHM for 
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Figure 5. Integrated Hi spectrum of NGC 300. The fluctuations 
near i>lsr = kms" 1 are due to Galactic foreground emission 
that was not deconvolved. The vertical dotted line indicates the 
systemic velocity of NGC 300 of i>lsr = 136 kms' 1 (equivalent 
to 144 kms -1 in the barycentric reference frame). 



Gaussian line profiles). These values are n ot atypical for the 
warm neutral medium of the Milky Way (|Kalberla fc Kerq 
2009) and indicate an upper limit on the order of 10 4 K 
for the gas temperature. The true gas temperature may be 
well below this upper limit, as effects such as turbulent mo- 
tion and beam smearing will contribute to the observed line 
width. 

Fig. [7] reveals regions of increased velocity dispersion 
near the boundary between the inner and outer Hi disc, 
in particular on the north-eastern side of the galaxy. The 
highest dispersions of about 20 kms -1 are found along a fil- 
ament in the eastern part of NGC 300. Such high values can 
no longer be explained by thermal line broadening alone, 
as in that case gas temperatures of about 5 x 10 4 K would 
be required. Instead, the high velocity dispersion suggests 
the presence of different gas components along the line of 
sight as supported by the spectra shown on the right-hand 
side of Fig. [7] The spectra were extracted from the loca- 
tions of highest dispersion marked with A and B in the map 
and reveal a clear double line profile at the boundary be- 
tween the inner and outer Hi disc. Apparently, each disc 
component produces a separate line component in the spec- 
trum, and the transition between the inner and outer disc 
is not gradual but abrupt within the size of our synthesized 
beam. This is also suggested by the radial velocity field as 
indicated by the white contours in Fig. [7| In the transition 
region between the inner and outer disc we observe a strong 
discontinuity in the velocity field. The resulting kinks in the 
velocity contours coincide precisely with the regions of in- 
creased velocity dispersion in the map, thus confirming our 
impression of an abrupt transition between the inner and 
outer disc of NGC 300. 

Interestingly, the transition from the inner disc to the 
outer disc is less sharp on the south-western side of NGC 300 
where the velocity contour lines are much smoother, and - 
despite a slightly increased velocity dispersion - the spectra 
do not exhibit a double-peak profile. Instead, there is a grad- 
ual transition from the inner to the outer disc component in 
the spectrum across several beam widths. 




19.0 - 
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Figure 6. H I column density profile of NGC 300 as a function of 
declination offset (solid curve). The uncertainties are plotted as 
the dotted curves and reflect our 5cr H I column density sensitivity 
of about 10 19 cm -2 . 

3.2 Rotation curve 

In order to determine the rotation curve, v Iot (r) , of NGC 300 
we must assume a model for the three-dimensional kinemat- 
ics of the galaxy. One of the simplest and most commonly 
applied models is to assume that the gas particles move in 
circular orbits along tilted ri ngs with varying position an- 
gle, u >, and inclination, i (e.g. iRogstad. Lockhart fc Wrightl 
Il974h . In addition, each ring is assigned a position, (xo, J/o), 
and systemic velocity, v sys , both of which can also be allowed 
to vary. 

We used the GIPSY task rotcur to fit a tilted ring 
model to the radial velocity field of NGC 300 as derived from 
the fitting of Gauss-Hermite polynomials to the spectral 
lines. We fitted 20 rings with radii ranging from r = 100 to 
2000 arcsec (0.9 to 18.4 kpc) in steps of 100 arcsec. Each 
ring was chosen to be 100 arcsec wide to match the angular 
resolution of our data. Data within an angle of ±20° around 
the minor axis of NGC 300 were excluded from the fit, and 
all data were weighted with the cosine of the position angle, 
|cob(0)|. 

We first derived the position of the centre of NGC 300 
from an optical l/-band image t aken with the MPG/ESO 
2.2-m telescope at La Silla l|Pietrzvhski et all 120021 ). 
NGC 300 has a very distinct, point-like core that is thought 
to be associated wi th either a nuclear stellar cluster or e ven a 
central black hole (|Soffner et al.lll99rj ; lKim et al.ll2004T i and 
should therefore provide us with an accurate position of the 
dynamical centre. A Gaussian fit to the optical core revealed 
a central position of a = 00 h 54 m 53!4 and S = -37°4l'02'.'6 
in the J2000.0 coordinate systemQ Next, we ran ROTCUR 
with the optically derived central position fixed and the gas 
expansion velocity also fixed to a value of f cxp = kms . 
All other parameters (systemic velocity, rotation velocity, 
position angle, and inclination angle) were still left as free 
parameters. After the first run of ROTCUR we determined the 
systemic velocity by averaging over the fitted systemic veloc- 
ities of all rings. The resulting value of v Bys = 136±2 kms -1 
in the LSR frame, corresponding to 144 ± 2 kms' 1 in the 

2 This position is consistent with the dynamical centre of the 
innermost tilted rings as derived by ROTCUR. 
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Figure 7. The greyscale image and black contour lines show the velocity dispersion of the H I gas in NGC 300 as derived from Gauss- 
Hermite polynomials fitted to the spectral lines. The contour levels are 10, 12, 14, 16, 18, and 20 kms -1 . The white contour lines show 
the velocity field for comparison and have been drawn in intervals of 20 kms -1 centred about the systemic velocity of dlsr = 136 kms -1 
(bold contour line; equivalent to 144 kms - in the barycentric reference frame). The two spectra were extracted from the regions of 
highest dispersion marked with the white circles labelled A and B in the map. 



barycentric reference frame, was then used and kept fixed 
in all subsequent runs. This result is consistent with the 
barycentric velo c ity o f v ByB = 146 ± 2 kms -1 found by 
iKoribalski et all l|2004h based on HIPASS. 

Next, we ran ROTCUR again with only the rota- 
tion velocity, position angle, and inclination angle left as 
free parameters. W e followed the approach described by 
Ide Blok et alJ lj2008T l and first fixed the radial position angle 
profile. To reduce the initial bumpiness of the position angle 
profile we smoothed it by applying a boxcar filter with a 
width of 300 arcsec (equal to three radial bins) to the posi- 
tion angle solution. After the next run of ROTCUR we also 
smoothed and fixed the inclination angle profile using the 
same boxcar filter as for the position angle. With all other 
parameters fixed, we then ran ROTCUR for the last time to 
obtain the final solution for the rotation curve. 

Our final Hi rotation curve of NGC 300 is shown in 
Fig. [gji and Table [2] A comparison between the observed 
and modelled velocity fields is presented in Fig. [9] Due to 
the combination of a large field of view with a fairly compact 
array configuration we are able to determine the rotation 
curve out to an angular radius of 2000 arcsec equivalent to a 
physical radius of 18.4 kpc. This is a significant improvement 
compared to the previous measurem ent with the Very Large 
Array (VLA) bv lPuche et all (|l990h which covered a radius 
of onl y 11 kpc. The rotation curve derived bv lPuche et al.l 
( 1990) is also plotted in Fig.[HJi for comparison. Within their 
uncertainties both rotation curves agree quite well, although 
a direct comparison is difficult due to the different resolu- 
tions of the VLA observations and the ATCA observations 
presented here. At a higher resolu t ion of 50 arcsec FWHM, 
the rotation curve of iPuche et ail l|l990h shows small-scale 
variations not present in our ATCA rotation curve. In ad- 
dition, there is some discrepancy in the very inner region of 
r < 5 arcmin which is most likely caused by different de- 
grees of beam smearing in combination with a rapid change 
in rotation velocity in the inner part of the galaxy. Hence, 



our ATCA observations in particular are not capable of reli- 
ably defining the rotation curve of NGC 300 across the inner 
3 kpc. 

The derived rotation curve of NGC 300 rises out to 
a radius of about 15 to 20 arcmin (equivalent to 8 to 
11 kpc) where we observe a maximum rotation velocity of 
Wrot = 98.8 ± 3.3 kms" 1 at r = 16.7 (9.2 kpc). T his value is 
consistent with the results of lPuche et a l. (1990) who found 
a maximum veloc ity of p rot = 97 . ± 5 .7 kms -1 at a radius 
of 15.'8 (8.7 kpc). lRogstad et all (|l979h derived similar val- 
ues in the range of iVot = 94 to 102 kms -1 at a radius of 
16 arcmin (8.8 kpc). Beyond about 10 kpc the rotation ve- 
locity gradually decreases to below 85 k ms" . This decrease 
was not seen in the previous VLA data o f lPuche et all £l990) 
due to insufficient radial extent of their map. 

A particular problem of rotation curve fits in general is 
the determination of uncertainties. The statistical errors of 
the fit, as provided by ROTCUR, are usually much smaller 
than the actual uncertainties of the solution which are dom- 
inated by systematic errors such as non-circular motion of 
the observed gas, insufficient spatial resolution of the ob- 
servations, effects of varying optical depth, etc. There are 
two possible ways to obtain more reliable estimates of the 
uncertainties. One possibility is to determine the standard 
deviation of velocities along each tilted ring from the mean 
value of the rotation velocity for that ring. This provides 
us with a more realistic assessment of the true deviations 
of the observed velocities from the model. The error bars 
of our final rotation curve presented in Fig. [8ji have been 
determined by this method. 

Another possibility is to derive the rotation curve for 
the approaching and receding side of the galaxy separately 
and assess the differences between the two curves. This is 
shown in Fig. [Hp where the approaching and receding curves 
are compared to the overall rotation curve. The fits for 
the approaching and receding sides are independent from 
the overall rotation curve fit, as we derived separate solu- 
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Figure 8. Panel a shows the derived overall rotation curve of NGC 300 out to a radius of 2000 arcscc. The rotation curve determined 



bv lPuche et al ] l|l990h based on VLA data is shown for comparison. Panel b shows the rotation curves derived from the approaching and 
receding side of NGC 300 separately compared to the overall rotation curve. Panels c and d show the position angle, <p(r), and inclination 
angle, i(r), respectively, of the tilted ring model used for the final rotation curve fit with ROTCUR (filled black circles) compared to the 
original values before smoothing (open grey circles with error bars). The dashed and dotted grey lines show the values derived separately 
for the approaching and receding side, respectively. 



tions for the position angle and inclination. The two rota- 
tion curves both show a clear maximum and then decrease 
again at larger radii. The position of the maximum as well 
as the overall shape of the rotation curve, however, differ 
substantially between the approaching and receding side. 
These discrepancies indicate systematic uncertainties of up 
to ±10 kms" 1 for the overall rotation curve of NGC 300. 
The rotation curve of the receding side continually decreases 
after reaching a maximum at r ~ 10 kpc, whereas the curve 
derived from the approaching side also decreases but then 
flattens out at larger radii. In addition, the approaching side 
has its rotation curve maximum significantly further inward 
at r ~ 6.5 kpc. These differences indicate a high degree of 
asymmetry in the kinematics of the gas disc of NGC 300, 
some of which is immediately evident from the velocity map 
in the left-hand panel of Fig. [§] The results indicate that 
the finding of a decreasing rotation curve in the outer disc 
should be approached with caution because non-circular mo- 
tions and asymmetries in the disc could have contributed 
significantly to the observed velocity field of NGC 300. 

Fig. fSJi and|H]i show the position angle, <p(r), and incli- 
nation, i(r), used in the final rotation curve fit with ROTCUR. 
For comparison we also show the originally fitted values be- 
fore applying the boxcar filter for smoothing and the sepa- 
rate solutions for the approaching and receding side of the 
galaxy. The position angle is constant over the entire inner 
disc of NGC 300 with <p(r) = 290° for r < 10 arcmin. Be- 
yond the edge of the inner disc the position angle of the 
gas disc gradually changes by about 40° only to become 



constant again in the very outer disc with <fi(r) = 332° for 
r > 30 arcmin. This behaviour of ip(r) reflects the visual 
impression of a twist in the Hi disc of NGC 300 and sug- 
gests the presence of two distinct gas discs, each with its 
own spatial orientation. The variation in inclination, i(r), 
with radius appears less regular. The inclination of the gas 
disc increases from about 40° to just over 50° within the 
inner 15 arcmin (8 kpc) and then gently decreases again to 
i(r) < 45° in the outermost regions of the disc. 

The observed decrease in rotation velocity is strongly 
coupled with the characteristics of the inclination angle in 
the outer parts of the disc. A smaller inclination angle would 
result in a higher rotation velocity. In order to obtain a flat 
rotation curve, the inclination angle of the outermost tilted 
ring would need to be about 35° . This is significantly smaller 
than the inclination angle of about 43° resulting from the 
tilted ring model. At the same time, we estimate a geometric 
inclination angle of about 44° by comparing the major and 
minor axis of the 2 x 10 19 cm -2 H I column density contour 
level under the assumption of perfect circular symmetry of 
the disc. This is consistent with the kinematic inclination 
from the tilted ring model. 

It is interesting to note in this context that the maxi- 
mum rotation velocity occurs near the transition radius be- 
tween the inner and outer H I disc. A closer inspection of 
the rotation curves derived for the approaching and reced- 
ing sides of NGC 300 (Fig. [8]) reveals that the position of 
maximum rotation velocity is closely linked to strong varia- 
tions in inclination angle at that particular radius. It is pos- 
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Figure 9. Panel a shows the observed radial velocity field of NGC 300 based on the position of the intensity maximum of Gauss- 
Hermite polynomials fitted to the spectral line profiles. Panel b shows the modelled radial velocity field derived from the overall rotation 
curve as shown in Fig. [8] In both panels, the contour lines are drawn in intervals of 15 kms — 1 centred about the systemic velocity of 
^LSR = 136 kms -1 (bold contour line; equivalent to 144 kms~ 1 in the barycentric reference frame). All velocities in the two maps are 
given in the LSR frame. Panel c shows the difference between the observed and modelled velocity fields. The contour lines are drawn at 
levels of ±4, ±8, ±12, and —16 kms -1 . The black, dotted ellipse in all three maps marks the outer edge of the model and illustrates the 
significant lopsidedness of the gas disc of NGC 300. 



sible that part of this variation in inclination is an artefact 
caused by the failure of the tilted ring model to cope with 
the sudden transition between the inner and outer disc. Con- 
sequently, the decrease in rotation velocity across the outer 
disc of NGC 300 could be less pronounced than indicated 
by the tilted ring model. 



4 MASS MODELS 

After having determined the overall rotation curve, we can 
now try to derive mass models of NGC 300 using the GIPSY 
task ROTMAS. This will allow us to study the radial distri- 
bution of the different visible and dark matter components 
in NGC 300 and assess how much dark matter is required to 
explain the rotation curve. For this purpose, we assume the 
presence of three mass components, namely a gaseous disc, 
a stellar disc, and a dark matter halo. ROTMAS will attempt 
to fit the velocity curves of these three components to the 
total rotation curve: 

"retW = /gasfgasM + /*V*(r) + UdmW- ( 2 ) 

The factors /* and / gas are mass scaling factors for the stel- 
lar and gaseous disc, respectively. ROTMAS determines the 
velocity curve, v(r), of a mass component from its mass sur- 
face density, E(r). In the following Sections 14.11 and 14.21 we 
will derive the required mass surface densities of the gaseous 
and stellar disc. In Section 231 we will introduce the different 
dark matter models used in the mass modelling of NGC 300. 
Finally, Section 14.41 will discuss the outcome of the mass 
modelling. 

4.1 Gaseous component 

Deriving the radial mass distribution of the gaseous disc 
in NGC 300 is relatively straightforward. We first used the 
GIPSY task ellint to derive the radial H I column den- 
sity profile from the column density map, using the tilted 
ring parameters obtained from the rotation curve fit with 



ROTCUR. The column density profile, Nui{r), can then be 
converted into gas mass surface density via 

S g as(r) = fmnNniir) cos(i) (3) 

where mH = 1.674 x 10 -27 kg is the mass of a hydrogen 
atom, and / is a mass correction factor. For our analysis we 
assume a value of / = 1.4 to account for primordial helium 
expected to be associated with the neutral hydrogen gas in 
the disc of NGC 300. The factor cos(i) is required to de- 
project all column densities to face-on values. The resulting 
mass surface density profile of the gas disc is shown in Fig. 1101 
and Tabled 

Throughout the inner part of the gas disc the surface 
density is basically constant with E gas « 7 M© pc -2 (includ- 
ing the correction for helium). Beyond a radius of about 
6 kpc the gas mass surface density begins to decline. The 
outer profile can be reasonably well described by either 
an exponential decline, E gas oc exp(— r/h), with a radial 
scale length of h — 3.44 ± 0.15 kpc, or by a power law, 
E gas oc (r/kpc) a , with a slope of a = -2.38 ± 0.07. 

This method, of course, assumes that the H I gas is opti- 
cally thin which might not be the case for the densest regions 
along the spiral arms of NGC 300. As a consequence, column 
densities, and hence gas mass surface densities, in these re- 
gions would be underestimated. In addition, we have made 
the assumption that the gas disc is infinitely thin, and any 
significant thickness of the disc would result in errors in the 
determination of H I column densities, depending also on the 
inclination angle of the disc. 



4.2 Stellar component 

Deriving the stellar mass distribution in NGC 300 is signif- 
icantly more complicated. In order to determine the radial 
mass profile of the stellar component of NGC 300 and its 
contribution to the derived rotation curve we need to find 
a way to convert the optical or near-infrared flux density 
profile, S\(r), into stellar mass surface density, E*(r), via 
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Table 2. Modelling parameters of NGC 300. The columns denote: galactoccntric radius, r, in arcsec and kpc; derived rotation velocity, 
v TO t] tilted ring position angle, tp, and inclination, i; stellar mass surface density, E*, as derived from the Spitzer IRAC 3.6 \xm, 4.5 um, 
and combined data; gas mass surface density, S gas . 
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the stellar mass-to- light ratio, Ta, such that 

E*(r)~T A S A (r). (4) 

The index A indicates that both the mass-to-light ratio as 
well as the flux density are wavelength-dependent and are 
usually determined within a specific photometric band of 
wavelength A. 

By using near-infrared instead of optical data one can 
largely overcome the problem of absorption by dust in the 
measurement of the radial flux density profile of the stellar 
component. NGC 300 has been observed with the Infrared 
Array Camera (IR AC) on board the Spitzer Space Telescope 
l|Helou et al.l 120041 ) . and for this work we have used IRAC 
data in the 3.6 and 4.5 u.m bands to determine the stellar 
surface brightness profile. The nominal waveleng ths of these 
two b ands are 3.550 and 4.493 um, respectively l|Fazio et al.l 
l2004h . IRAC pipeline images are usually calibrated in terms 
of flux density in units of MJy sr~ . 

The flux density profile, S\(r), in the 3.6 and 4.5 u.m 
wavebands must then be converted to stellar mass density, 
£*(r), using the stellar mass-to-light ratio in units of the 
solar value, = T A /Y®. The corresponding conversion 
formula in Eq. [4] can be rewritten in the form 

E*(r) = C A T' A ^|M (5) 

•30 

where C\ is a wavelength-dependent conversion factor, /a 
is the wavelength-dependent aperture correction factor for 
extended sources in IRAC images, and So denotes the 
wavelength-de pendent zero- m agnitude flux. Based on the 
calculations of lOh et al.l |2008l ) we derive conversion factors 
in the 3.6 and 4.5 u.m wavebands of 

C 3 .6 =0.196 Mope" 2 , (6) 

C4.5 = 0.201 Mope -2 . (7) 



According to iReach et al.l |2005h , the aperture correction 
factors for the 3.6 and 4.5 urn bands are /a = 0.944 and 
0.937, respectively, and the corresponding zero-magnitude 
fluxes are So = 280.9 ± 4.1 and 179.7 ± 2.6 Jy. 

Before we can use Eq. \S\ to convert radial flux density 
profiles into stellar mass density profiles we first need to 
determine the stellar mass-to-light ratio in NGC 300. De- 
termination of the stellar mass-to-light ratio is very difficult 
because it depends on several uncertain parameters such as 
the initial mass function (IMF) of NGC 300, its star forma- 
tion history, metallicity, etc. To make matters worse, most 
of these parameters will change with growing distance from 
the centre of a galaxy, resulting in a radial variation, T\(r), 
of the stellar mass-to-light ratio. 

Colour indices provide a convenient way to determine 
the m ass-to-light ratio throu gh basic photometric measure- 
ments. iBell fe de Jond (|200lh used galaxy evolution models 
to explore the relation between stellar mass-to-light ratio 
and optical/near- infrared colours of galaxies. They found a 
strong and robust correlation between the stellar mass-to- 
light ratio and the integrated colour of the stellar population 
in t heir model gala xies. They also concluded that a modi- 
fied |SaIp_eterJ (1 19551 ) IMF with a smaller fraction of low-mass 
stars is consistent with the observed stellar mass-to-light ra- 
tios of a sample of spiral galaxies used for comparison with 
the models. From the relations between ma ss-to-light ratio 
and colour found by IBell fe de Jond l|200lf ) we derive the 
following relation between the A'-band stellar mass-to-light 
ratio, T' K , and the J — K colour index: 

log 10 (r' K ) = 1.434 (J-K)- 1.380. (8) 

The relation was derived using the mass-dependent forma- 
tion epoch model with bursts and the modified Salpeter IMF 
adopted bv lBell fc de Jond |200ll ). Since our surface bright- 
ness measurements are based on Spitzer data, we have to 
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Figure 10. Azimuthally averaged radial profiles of the stellar 
mass surface density, S*(r), and the gas mass surface density, 
S gas (r), of NGC 300 (see Table 0. 



convert T' K into the stellar mass-to-light ratios derived for 
the Spitzer IRAC bands of 3.6 and 4 . 5 urn . Using stellar 
population synthesis models, foh et al.l (2008) found the fol- 
lowing well-defined correlations: 



Tg.g = 0.92 T' K - 0.05, 
T^.s = 0.91 Tk -0.08. 



(9) 
(10) 



With the set of equations derived above we now have all 
tools at hand to determine the stellar mass surface density 
of NGC 300. 

From the Two Micron All Sk y Survey (2MASS) Large 
Galaxy Atlas jjarrett et alj|2003 ) we derive a global colour 
index of J - K = 0.646 ± 0.033 mag for NGC 300. Inserting 
this value into Eq. [8] yields a stellar mass-to- light ratio for 
NGC 300 in the 7f-band of 



log 10 (T' K ) = -0.454 ± 0.047. 



(11) 



Using Eg. l9l and 1 101 we derive stellar mass-to-light ratios in 
the IRAC 3.6 and 4.5 um bands of 



log 10 (T!,. 6 ) = -0.564 ±0.055, 
logi (T4.5) = -0.620 ± 0.054. 



(12) 
(13) 



We can now use Eq. [5]to derive the stellar mass surface den- 
sity of NGC 300. Again, we used the GIPSY task ellint 
with the tilted ring fitting parameters derived earlier to de- 
termine the surface brightness profiles of NGC 300 from the 
Spitzer 3.6 and 4.5 pm images. Both profiles were then cor- 
rected for inclination and combined into a single stellar mass 
surface density profile which is shown in Fig.[l0]and Table(2] 
The resulting total stellar mass is (1.0 ±0.1) x 10 9 M . 

The derived profile can be described by an exponen- 
tial decline, £*(r) oc exp(— r/h), with a radial scale length 
of h = 1.39 ± 0.04 kpc. This value is in good agreement 
with the scale len gth of 1.38 kpc measured in the / band by 
iKim et all (|2004T ). Comparison w ith the values of 2.06 kpc 
in the Bj band (ICarignanl 119851 ) and 1.94 kpc in the B 
band l|Kim et "all I2004T) suggests that the outer stellar disc 
of NGC 300 is bluer than the inner disc. 

As mentioned earlier, we cannot expect the stellar mass- 
to-light ratio to be constant across the entire galaxy, but 
there will likely be a radial variation in Ta reflecting differ- 



ent star formation histories and stellar populations at differ- 
ent distances from the centre of NGC 300. A way to account 
for this radial variability of Ta would be to determine the 
radial profile of the J — K colour index of NGC 300 based 
on infrared imaging data. Unfortunately, the signal-to-noise 
ratio in the 2MASS images is very low, and determination 
of Ta as a function of galactocentric distance is therefore 
not feasible due to large uncertainties. 



4.3 Dark matter component 

4-3.1 Pseudo-isothermal halo 

The spherical pseudo-isothermal halo is one of the simplest 
models for the density profile of dark matter haloes. Its den- 
sity profile resembles that of an isothermal sphere and reads 

Qo 



g(r) = 



(14) 



1 + (r/r c ) 2 

(e.g. iBegeman. Broeils fc Sanderj|l99ll ) where go is the cen- 
tral density of the halo and r c the so-called core radius. The 
corresponding velocity profile is 



v 2 (r) = 4nGg r 2 
(|Kentlll986l 1. 



r c IT 

1 arctan — 

r \r c 



(15) 



4.3.2 Burkert halo 

For his study of dwarf galaxy rotation curves iBurkertl (|l995f) 
introduced the following empirical halo density profile: 

„3 



g(r) = 



gor c 



(r + r c )(r 2 ± r 2 ) ' 



(16) 



Again, go and r c are the central density and core radius of 
the halo, respectively. The resulting velocity profile reads 

6.4G£>or c 3 



v 2 (r) 



In 1 + 



1 + 



arctan 



(17) 



l|Salucci fc Burkertft OQO). The major difference between the 
Burkert halo and the pseudo-isothermal halo is that for 
large radii the Burkert halo density behaves as g(r) ~ r -3 , 
whereas for the pseudo-isothermal halo we get a different 
behaviour of g(r) ~ r~ 2 . 

4.3.3 NFW halo 

The NFW dark matter density profi le was suggested by 
iNavarro. Frenk fc White] j 19951 . 1996) based on numerical 
simulations of dark matter haloes in a hierarchically clus- 
tering universe. They found that in their models the radial 
density profile of dark matter haloes can be accurately de- 
scribed by 

5c£>crit 



g(r) = 



7- 1 + 



(18) 



where g CT \t is the critical density of the universe, S c is the so- 
called characteristic, dimensionles s density of the halo, and 
r s is a scale radius. According to iNavarroT Frenk fc White] 
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Table 3. Results of the mass modelling of NGC 300 for the pseudo-isothermal (ISO), Burkert, and NFW dark matter halo profiles. /+ 
and /gas are the mass scaling factors for the stellar and gaseous disc, respectively. The parameters in columns 4 and 5 are the free model 
parameters of the different dark matter models as explained in the text. Column 6 lists the reduced \ 2 value of the model as a measure 
of the goodness of fit. The last two columns list the total mass, Mtot, of the model within a radius of r = 18.4 kpc and the fraction of 
dark matter, /dm- Note that the NFW fit with fixed gas mass produced negative stellar masses and is therefore not listed. 



Halo model 


/* 


/gas 




eo 


-^red 




Mtot 


/dm 








(kpc) 


(M oP c- 3 ) 




(10 10 


M Q ) 




ISO (fixed) 


1 


1 


0.93 ±0.14 


0.170 ±0.046 


2.10 




3.6 


0.92 


ISO (gas fixed) 


1.6 ±1.3 


1 


1.17 ±0.61 


0.107 ±0.111 


2.16 




3.6 


0.90 


ISO (free) 


2.3 ±0.6 


4.8 ± 1.0 


0.74 ±0.35 


0.135 ±0.119 


1.17 




3.3 


0.65 




/* 


/gas 




eo 


Xred 




Mtot 


/dm 








(kpc) 


(M 0P c- 3 ) 




(10 10 


M Q ) 




Burkert (fixed) 


1 


1 


2.78 ±0.14 


0.081 ± 0.008 


1.18 




3.2 


0.91 


Burkert (gas fixed) 


2.2 ±0.9 


1 


3.58 ± 0.75 


0.044 ±0.021 


1.20 




3.3 


0.88 


Burkert (free) 


3.1 ±0.9 


3.0 ± 1.4 


4.04 ± 1.10 


0.025 ±0.018 


1.12 




3.3 


0.73 




/* 


/gas 




T200 


^red 




M tot 


/dm 








(kpc) 


(kpc) 




(10 10 


M ) 




NFW (fixed) 


1 


1 


5.81 ± 0.56 


89.7 ±2.4 


1.47 




3.4 


0.91 


NFW (free) 


1.1 ± 1.1 


3.9 ± 1.0 


4.00 ± 1.53 


71.9 ±6.7 


0.95 




3.3 


0.74 



l| 19971 ). the resulting circular velocity curve of the dark mat- 
ter halo reads 



V(r> * ln(l + c)-lfs 



(19) 



with the substitution x — r/r2oo and the concentration pa- 
rameter c = r-2oo/fs (not to be confused with the speed of 
light, c). In this form of the equation, r2oo denotes the radius 
within which the mean density of the halo equals 200 times 
the critical density of the universe, g C iit = 3Hg/(87tG). This 
choice is somewhat arbitrary, although r2oo corresponds to 
approximately the virial radius of the halo. Accordingly, V2pp 
is the circular velocity of the halo at V2oo expressed by 



/ GM200 1 
^200 = \l = 10Hor2oo. 

»"200 



(20) 



This leaves us with only two free parameters for the rota- 
tion curve fit, namely r a and T2oo- In additio n, we assumed a 
Hubb le constant of Ho = 70.5 kms -1 Mpc -1 (Hinshaw et al. 
120091 ) throughout this paper. Note that the NFW density 
profile has a singularity at the centre of the halo. Neverthe- 
less, both the mass and circular velocity are well-defined and 
finite at all radii. 



4.4 Mass modelling of NGC 300 

Throughout the literature there has been no consistent ap- 
proach as to what vertical de nsity profile, g(z), to assume for 
the ga seous and stellar discs. ICapozziello. Cardone fc Troisil 
l|2007l) . for ex ample, assumed bot h discs to be infinitely 
thin, whereas Ide Blok et al.l ()2008h assumed an infinitely 
thin gas disc but a sech 2 (z/zp) dist ribution for the stellar 
disc. ISanders fc Noorderm ecr (2007), in turn, also assumed 
the gas disc to be extended in the vertical direction with the 
same thickness as that of the ste ll ar dis c. In contrast to all 
these approaches, iKalberla et all (|2007m demonstrated that 
the H 1 disc of the Milky Way shows strong flaring, with the 
scale height (half width at half maximum) increasing from 



zp = 60 pc at a radius of 4 kpc up to zo w 2.7 kpc at a 
radius of 40 kpc. 

The sech 2 (z/zo) distribution of the stellar disc is 
motivated by studies of e dge-on spiral galaxies by 
van der Kruit fc Searld (|l981al lbf). They found that the mass 
density of the stellar disc is of the form 



g(r,z) 



£(r)sech (z/zo) 



(21) 



with the s cale height, zp, being l argely independent of radius, 
r (but see Ide Griis fc Peletierl l|l997h for evidence of a vari- 
ation of zo with r). We adopt this distribution in our model 
of NGC 300 and assume a ratio of h/zp = 5 in accordanc e 
with previous studies (e.g. Ivan der Kruit fc Searld Il981ah . 
From the exponential scale length, h — 1.39 ± 0.04 kpc, of 
the combined 3.6 and 4.5 urn mass surface density profile 
we therefore derive a vertical scale height of zo ~ 280 pc. 

As for the gaseous disc, we have considered two cases, 
namely an infinitely thin disc and an exponential distribu- 
tion with a scale height equivalent to that of the stellar disc. 
It turns out that a thicker gas disc results in a slightly lower 
rotation velocity and a general smoothing of the velocity 
profile of the gaseous component. However, the introduc- 
tion of a non-zero scale height for the gaseous disc does not 
have any significant influence on the mass modelling as long 
as the scale height is not larger than a few hundred pc. We 
therefore decided to make the assumption of an infinitely 
thin gas disc in all our models. 

For all three halo models we carried out fits for fixed, 
partly fixed (gas only) , and variable stellar and gaseous mass 
components! 3 ] The resulting mass model fits are shown in Ta- 
ble [3] and Fig. 1111 All three dark matter halo models result 
in reasonably good fits. The goodness of fit of the pseudo- 
isothermal halo model is somewhat lower compared to the 



3 Note that the fit for the NFW halo with fixed gas mass scaling 
factor failed, because the best-fitting model resulted in a nega- 
tive stellar mass scaling factor. Therefore, this case will not be 
presented and discussed here. 
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Table 4. Results of the mass modelling of NGC 300 for New- 
tonian dynamics and Modified Newtonian Dynamics (MOND) 
under the assumption of no additional dark matter halo. /* and 
/gas are the mass scaling factors for the stellar and gaseous disc, 
respectively. The last column lists the reduced \ 2 value of the 
model as a measure for the goodness of fit. 



Dynamics 




/* 


/gas 


^rcd 


Newtonian (gas fixed) 


9.5 ± 


1.4 


1 


67.42 


Newtonian (free) 


4.2 ± 


0.3 


8.4 ±0.3 


2.17 


MOND (gas fixed) 


2.2 ± 


0.2 


1 


3.31 


MOND (free) 


2.6 ± 


0.2 


0.6 ±0.2 


2.50 



Burkert and NFW halo models. The lower performance of 
the pseudo-isothermal halo is mainly due to its r~ 2 den- 
sity profile at large radii. As a consequence, the velocity 
curve will continue to increase and eventually converge to a 
parameter-dependent constant, Voo- Therefore, the pseudo- 
isothermal halo model cannot explain the decreasing rota- 
tion curve of NGC 300, whereas the Burkert and NFW halo 
models (both with an r -3 behaviour) are doing better. Of 
course, all three halo models produce a lower reduced \ 2 in 
the case of free stellar and gas mass scaling factors. Whereas 
for the pseudo-isothermal and NFW haloes the improvement 
is significant, the Burkert halo is the only model to produce 
a good fit with fixed mass scaling factors. In the case of 
free mass scaling factors we obtain values in the range of 



/* w 1 ... 3 and / g; 



. 5. Fixing the gas mass scal- 



ing factor to / gas = 1 yields somewhat smaller stellar mass 
scaling factors of /* w 1.6 ... 2.2. 

As the stellar mass scaling factor strongly depends on 
the assumptions made about the stellar mass-to-light ratio 
of NGC 300, a higher value of /* can arise from differences in 
star formation history or metallicit y compared to the stan - 
dard spiral galaxy model found by iBell fc de Jond l|200ll ) . 
In comparison, the best mass model fit for NGC 300 found 
by iKenti l|l987l ) yields /* = 1.0 (relative to our stellar disc 
mass of 1.0 x 10 9 Mm), and his maximum-disc model has 
/* = 3.1. iPuche et alj(|l990f i derive somewhat larger values 



of /* = 2.2 and 4.3 for their best-fitting and maximum- 
disc models, respectively. These numbers are all in general 
agreement, but they also expose the uncertainties in the de- 
termination of the stellar mass-to-light ratio and the total 
stellar disc mass of NGC 300. 

The free mass scaling factors for the gaseous component 
are somewhat higher than expected. A factor of 2 can still 
be explained by optical depth effects in the inner regions of 
the gas disc, significant amounts of ionised gas throughout 
the disc, or missing flux in our Hi map due to lack of short 
interferometer spacings. Gas mass scaling factors of 4 or 5, 
however, seem unreasonably large and not realistic. 

Note that in the very inner region of NGC 300 the 
value of the rotation velocity of the gas disc is negative 
(see Fig. 1111) . This does not imply counter-rotation of the 
gas or the presence of negative mass, but simply reflects a 
net outward force on test particles due to the central de- 
pression of the gas disc. This net outward force results in a 
negative value of v 2 ^ and therefore an imaginary velocity in 
the mathematical description of mass models that is usually 
expressed by negative values of w gas . 



In addition to the dark matter models discussed above 
we also fitted two models without any dark matter compo- 
nent at all, using Newtonian dynamics a nd (non-relativ istic) 
Modified Newtonian Dynamics (MOND; [Milgrom|[l983T ). re- 
spectively. In both cases, only the stellar and gas mass scal- 
ing factors were left as free parameters of the fitting proce- 
dure. For the MOND model we assumed a constant qp = 
1.21 x 10" 10 ms" 2 l|Begeman. Broeils fc Sanders! ll99ll )FI 
The resulting mass models are shown in Tableland Fig. 1121 

Neither of the two models produces meaningful results. 
The best fit in the MOND scenario (\ 2 cd = 2.50) yields 



mass scaling factors of /* = 2.6 ± 0.2 and / g; 



0.6 ±0.2 



for the stellar and gaseous disc, respectively. Apart from the 
significantly larger Xred compared to all of the dark mat- 
ter models with free mass scaling factors, the MOND model 
results in a considerable downscaling of the gas disc mass. 
This is an implausible result because the derived gas mass 
surface density was directly derived from the observed H I 
column densities and must therefore be a lower limit. Fix- 
ing the mass scaling factor for the gas disc to a value of 
/gas = 1 slightly reduces the stellar mass scaling factor to 
/* = 2.2 ± 0.2, the goodness-of-fit parameter, however, sig- 
nificantly deteriorates to Xred = 3.31. 

The reason for the inferior performance of MOND is 
the decrease in rotation velocity at larger radii. MOND has 
been designed to explain the flat rotation curves originally 
observed in many galaxies. Therefore, the model will fail 
in cases of decreasing rotation curve. Consequently, deep H I 
observations with the aim to detect the faint, outer gas discs 
of galaxies have the potential to challenge MOND through 
the detection of decreasing rotation curves in the outer re- 
gions of many more galaxies. It should be noted, though, 
that - as discussed earlier - the observed decrease in the 
rotation velocity of NGC 300 is uncertain due to the possi- 
bility of non-circular motions and asymmetries in the disc. 
Deep studies of a much larger sample of spiral galaxies will 
be required to investigate the frequency of occurrence of de- 
creasing rotation curves in these galaxies. 

The Newtonian scenario without dark matter also fails 
to produce meaningful results. The best fit again has a large 
X 2 c d = 2.17. In addition, both the stellar and gaseous disc 
masses would have to be scaled up significantly by scaling 
factors of /* = 4.2 ± 0.3 and / gas = 8.4 ± 0.3, respectively. 
Such large scaling factors, in particular for the gas disc, are 
highly implausible. Of course, fixing the gas mass scaling 
factor to /gas = 1 does not produce any meaningful result at 
all, because the remaining free stellar mass component dra- 
matically fails to reproduce the overall shape of the rotation 
curve, resulting in a Xred of 67.4. 

In summary, we obtain the best mass model fits under 
the assumption of a Burkert or NFW dark matter halo. Both 
dark matter models can cope with the decreasing rotation 
curve of NGC 300 at larger radii. The Burkert halo provides 
the best fit with fixed mass scaling factors, but both halo 
models produce similarly good fits with variable mass scaling 
factors, although the gas mass scaling factor is slightly larger 
and less realistic for the NFW halo. The resulting total mass 
of NGC 300 within 18.4 kpc amounts to (2.9±0.2) x 10 10 M Q . 



4 Fits with variable ao and fixed mass scaling factors resulted in 
very similar values of ao ~ 1.3 X 10" 
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Figure 11. Results of the fitting of mass models to the rotation curve of NGC 300. The panels show, from top to bottom, the results 
for the pseudo-isothermal halo, the Burkert halo, and the NFW halo. For each halo model the results for fixed, partly fixed (gas only), 
and free stellar and gaseous mass components is shown. The black data points with error bars show the observed rotation curve. The 
stellar, gaseous, and dark matter components are plotted as the dotted, dashed, and solid black lines, respectively. The total velocity 
profile of the model is plotted as the solid grey line. Note that the NFW fit with fixed gas mass produced negative stellar masses and is 
therefore not displayed. 



For fixed mass scaling factors almost 90 per cent of the mass 
is contributed by dark matter. In the case of variable mass 
scaling factors we obtain 68 per cent dark matter for both 
the Burkert and NFW halo models (but only 57 per cent for 
the pseudo-isothermal halo model). 

However, all these results should be approached with 
great caution. First of all, none of the dark matter models 
discussed above is physically motivated. Instead, these are 
heuristic models based mainly on the results of numerical 
computer simulations. Secondly, the lower Xred values of the 
models with dark matter compared to those without dark 
matter may simply reflect the larger number of free fitting 
parameters in the former case. Any increase in the number of 
mass components or free parameters may naturally improve 
the goodness of fit, no matter what dark matter model one 
assumes. 



5 DISCUSSION 

5.1 Origin of the outer disc 

A conspicuous feature of NGC 300 is its extended outer H I 
disc spanning more than 1° (equivalent to about 35 kpc) 
across the sky. Although the outer disc was partly mapped 



by previous observations (e.g., IPuche et al.l Il990l ). this is 
the first time that its extent has been imaged out to the 
10 19 cm -2 column density level with moderately high spa- 
tial resolution of approximately 1 kpc. The results of our 
tilted ring model suggest that there is a substantial change 
in position angle between the inner and outer disc, resulting 
in a twisted appearance of NGC 300 which is particularly 
obvious in the velocity field. 

One possible scenario for this twist is that the distor- 
tion and warping of the outer disc of NGC 300 was caused by 
tidal forces during a recent encounter with another galaxy. A 
potential candidate for a close encounter could be NGC 55. 
Both galaxies are separated by about 8° on the sky, cor- 
responding to a projected separation of about 270 kpcQ 
Under the assumption that both galaxies move away from 
each other at a relative velocity of 200 kms -1 , there would 
have been a close encounter between NGC 55 and NGC 300 
about 1.3 Ga ago. This time scale is comparable to the 
1.5 Ga that have passed since the previous perigalactic pas- 
sage of the Magellanic Clouds on their orbit about the Milky 
Way and the resulting creation of the Magellanic Stream 

5 NGC 55 has ap proximately the same d istance from the Milky 
Way as NGC 300 JPietrzvriski et al.ll20od) . 
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Figure 12. Same as in Fig. 1111 but for non-relativistic MOND (top) and Newtonian dynamics (bottom). In both cases, no additional 
dark matter was included. In the left-hand column the gaseous mass component was held fixed whereas the stellar component was allowed 
to vary. In the right-hand column both components were variable. 



as predicted by numerical simulat ions ((Gardiner fc Noguchil 



2009). 



1996; Yoshizaw a fc Noguchill2003l ; however see iMastropietrol 



Alternatively, the distortion of the H I disc of NGC 300 
could have been caused by a recent encounter with a smaller 
companion galaxy. The nearest currently known companion 
is the dwarf-spheroi dal galaxy ESQ 294— G0 10. Its distance 
of about 1.9 Mpc l|Karachentsev et al.l [2003) is the same as 
that of NGC 55 and NGC 300, and its angular separation 
from NGC 300 is 6? 9, equiva lent to a physical sep aration of 
about 230 kpc. Its small size (|Bouchard et al.ll2005l ) and rela- 
tively large separation from NGC 300, however, suggest that 
ESO 294— G010 has never had any significant gravitational 
influence on NGC 300. At the same time, ESO 294-G010 is 
much closer to NGC 55 with an angular separation of only 
3? 5 equivalent to about 115 kpc. 

Deep optical observations of two fields in the outer 
part of NGC 300 with the G emini South 8-m telescope by 
iBland-Hawthorn et all (2005) revealed an extended stellar 
disc reaching out to a radius of at least 24 arcmin (equivalent 
to about 14 kpc or 10 optical scale lengths). Surprisingly, 
the luminosity of the stellar disc of NGC 300 shows a sim- 
ple exponential profile over its entire radial extent (except 
for the nuclear region) without any indication of a devia- 
tion or break. This seems to contradict our observation of a 
strongly warped H I disc with a rather discontinuous transi- 
tion between the inner and outer disc. If this warp were the 
result of tidal interaction, the stellar disc should have been 
affected in the same way. 

It is important to note that IBland-Hawthorn et al.l 

(2005) did not map the entire two-dimensional distribution 
of the outer disc on the sky, but observed only two fields in 
the south-eastern quadrant of NGC 300. Any warp in the 
stellar disc would therefore have to affect the stars at this 
particular azimuthal angle to be discernible, and further op- 



tical imaging at different azimuthal angles would be required 
to firmly exclude the presence of a warp in the outer stellar 
disc. 

IBland-Hawthorn et alj l|2005t ) discuss some of the impli- 
cations of this discrepancy and speculate about the origin 
of the extended stellar disc, including the possibility of stars 
being scattered from the inner disc. Another possibility is 
that gas accreted by NGC 300 through hot or cold accretion 
modes could have fuelled low levels of star formation in the 
outer disc. In this case, however, there should be a break in 
the radial stellar surface brightness profile. 

Alternatively, stars could have formed in the outer H I 
disc as a result of instabilities such as spiral density waves. In 
order to test this hypothesis, IBland-Hawthorn et al.l (|2005T l 
calculated the Toomre Q parameter, 



Q(r) = 



er(r)/c(r) 



(22) 



3.36GE gas (r) 

l|Toomrj|l96i ), based on the H I data of lPuche et afl l |l990j) . 

Here, a(r) is the velocity dispersion of the gas disc, n(r) is 
the epicyclic frequency of the rotating disc defined by 



r dr 



(23) 



and Egas(r) denotes the gas mass surface density, all 
three of which are a f unctio n of galactocentric radius, r. 
IBland-Hawthorn et all (2005) found values in excess of Q ~ 
5 for the disc beyond a radius of 10 arcmin (equivalent to 
about 5.5 kpc), suggesting that the outer disc of NGC 300 
is stable with respect to axisymmetric instabilities. 

Recalculation of Q based on our Hi observations with 
the ATCA, however, suggests somewhat smaller Q param- 
eters over much of the disc of NGC 300. Using Eq. 1221 and 
assuming t he same constant gas velocity dispersion of a — 
5 kms -1 as lBland-Hawthorn et all (|2005l ) yields typical val- 
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Figure 13. Toomre Q parameter of the gas disc of NGC 300 as 
a function of galactocentric radius. 



ues of Q as 1 ... 3 between 2 and 12 kpc radius fFig. H3[) , sug- 
gesting that the gas disc may be susceptible to axisymmet- 
ric instabilities. However, a slightly higher and more realis - 
tic velocity dispersion of 10 kms -1 (|Tamburro et al.ll2009l ) 
would result in the disc being Toomre stable at most radii. 
Furthermore, recent studies suggest that the star formation 
efficiency in disc galaxies is regulated by other effects, such 
as the physics of the ISM on sm all scales (jLerov et alj|2008h 
or the gas pressure in the disc l)Blitz fc Rosolowskvl f2006). 
and that the Toomre Q parameter does not necessarily pro- 
vide an accurate description of the ability of the disc to form 
stars. 

Obviously, the next logical step in our attempt to un- 
derstand the evolutionary history of NGC 300 would be a 
more detailed study of the outer stellar disc to determine 
the age and chemical composition of the stars. This was 
recent ly achieved by IVlaiic. Bland-Hawthorn fc Freeman! 
(2009) through g and i' photometry of three fields in the 
outskirts of NGC 300 with the Gemini South telescope. The 
observations confirmed that the stellar population in the 
outer disc is predominantly old. At the same time, there 
appears to be a radial gradient in metallicity with a sud- 
den break at a radius of about 10 kpc. While the inner disc 
shows a radially decreasing metallicity, this trend surpris- 
ingly reverses in the outer disc. Interestingly, the break oc- 
curs near the boundary between the inner and outer H I disc 
of NGC 300, although it is not clear at this point whether 
this is just a coincidence or whether there is a physical con- 

nect ion between the two phenomena. 

IVlaiic. Bland-Hawthorn fc Freeman! (|2009h propose 
two possible scenarios for the observed break in the metal- 
licity gradient. Either there could have been some kind of 
radial mixing of stars in the disc, or star formation could 
have slowly progressed outward during the evolution of 
NGC 300, resulting in a moderate metallicity gradient in 
the inner disc, but an almost flat and pristine metallicity 
distribution across the outer disc. In both cases, the gas 
in the outer disc would not play any significant role in 
the recent star formation activity of NGC 300, which is 
consistent with the generally old age (> 1 Ga) of the stellar 
population in the outer disc and the presence of only a 
small fraction of younger (a few hundred million years) 
main sequence stars. 



5.2 Ram pressure interaction 

The aforementioned strong asymmetries found in the outer 
Hi disc of NGC 300 are suggestive of ram pressure inter- 
action with the surrounding intergalactic medium (IGM) 
of the Sculptor Group. As NGC 300 is moving through 
the IGM, gas at the leading edge of the disc would be 
compressed, whereas gas at the trailing edge of the disc 
would be dragged and stripped by ram pressure forces 
(see, e.g., IVollmerl |2009| and references therein). The in- 
fluence of ram pressure has been mainly observed in the 
dense intergalactic medium of massive galaxy clusters (e.g. 
IWhite et all Il99ll ; IVollmer et "all 120041). in pa rticular the 
nearby Virgo cluster ( Giovanelli fc Havnesl 19831 ) , where ram 
pressure is believed to be res ponsible for the strip ping and 
removal of gas from galaxies (|Vollmer et al.ll200ll ). In con- 
trast, only very few cases of ram press ure stripping in galax 



alaxy 
2005; 



groups have been reported so far (e.g.,|K antharia et al, 
iMcConnachie et al. I l2007l ; iBouchard et al.|[2007l ). 

To assess whether ram pressure could have caused the 
observed asymmetries in the gas disc of NGC 300 we can 
compare the expected ram pressure forces with the gravita- 
tional forces stabilising the disc. The pressure imposed on 
the gas disc by a galaxy's own gravitational potential can 
be expressed as 



t)z 



(24) 



jGunn fc Gottl Il972l ; iRoediger fc Henslerl l2005h . Here, 
Sgas (r) is the gas mass surface density of the disc as a func- 
tion of radius, and $(r) denotes the gravitational potential 
of the galaxy. The term on the right-hand side of Eq. [24] 
determines the maximum of the absolute value of the gravi- 
tational acceleration perpendicular to the disc, with z being 
the height above the disc. 

The gas disc will be subject to ram pressure stripping 
if the pressure exerted by the IGM exceeds the stabilising 
gravitational pressure as determined by Eq. 1241 The ram 
pressure of the IGM is determined by the density of the IGM, 
£>igm, and the relative velocity of the galaxy with respect to 
the surrounding medium, « re i, thus 

-Pram = £>IGMt>rel- (25) 

Ram-pressure stripping of gas in the disc will occur when- 
ever the instability condition P ram > P gra v is fulfilled. As 
the gravitational pressure will be a strong function of galac- 
tocentric radius, one can also define a truncation radius at 
which ram pressure forces start to dominate over gravity. 
The simple analytic approach described above assumes that 
the galaxy is moving face-on through the IGM, however, hy- 
drodynamical simulations have demonstrated that even for 
larger inclination angles of the disc of up to about 60° with 
respect to the direction of motion the analytic approa ch pro- 
duces meaningful results (|Roediger fc Briigg en 2006). 

In order to assess the stability of the gas disc of 
NGC 300 against ram pressure stripping we have calculated 
the instability condition for different scenarios. We assumed 
the gravitational potential of the galaxy to be dominated by 
its dark matter halo and neglected the potential of the stel- 
lar disc and gas disc. This assumption will be met at greater 
heights, z, above the disc plane in the outermost regions of 
the disc where we expect the gas to be particularly vulner- 
able to ram pressure stripping. Furthermore, we assumed 
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Figure 14. Ratio between gravitational pressure and ram pressure in the disc of NGC 300 as a function of radius. The left-hand panel 
shows the pressure ratio for a constant velocity of v rei = 200 kms -1 and five different densities of (from top to bottom) njGM = 5 X 10 — 6 , 
1 x 10~ 5 , 2 x 10" 5 , 5 x 10" 5 , and 1 x 10~ 4 cm . The right-hand panel shows the pressure ratio for a constant density of rtiQM = 
5 x 10 -5 cm~ 3 and five different velocities of (from top to bottom) v Ie \ = 100, 150, 200, 250, and 300 kms -1 . The region above the 
horizontal black line at log 1 Q(P grav /P ram ) = is dominated by gravity, and the gas disc is stable. The region below that line is dominated 
by ram pressure, and the gas will be subject to ram pressure stripping. 



the dark matter halo to have an NFW profile with the pa- 
rameters (for the model with fixed stellar and gaseous mass 
scaling factors) derived in Section T4.4I and listed in Table [3] 
The gravitational potential of the NFW halo reads 

GMaoo Infl + f ) 
SnfwW = p i "Mr (26) 

r[ln(l + <0- (T+7)J 

l|Havashi. Navarro fc Spri ngel 2007) where M200, r s , and c 
are parameters of the NFW profile as discussed in Sec- 
tion [4X3] Solving Eq. [24] for the NFW potential and com- 
bining it with Eq. [25] allows us to investigate the instability 
condition for different values of the IGM density and relative 
velocity of NGC 300. The results are presented in Fig. 1141 

The left-hand panel of Fig. [TJ] shows the situation for 
a constant velocity of v TC \ — 200 kms -1 of NGC 300 with 
respect to the IGM but for five different IGM densities. The 
right-hand panel of Fig. [T3] shows the situation for a con- 
stant IGM density of 71igm = 5 x 10 -5 cm -3 but for five 
different velocities. As we would expect, most of the gas 
disc of NGC 300 is many orders of magnitude above the 
instability condition of log 10 (P grav /P ram ) < represented 
by the horizontal black line. However, under certain condi- 
tions the outermost part of the gas disc beyond a radius 
of about 15 kpc will get near or just below the instability 
line, in particular for IGM densities in excess of about a few 
times 10~ 5 cm -3 and sufficiently high relative velocities of 
Vrel 200 kms -1 . The density values assumed here are con- 
sistent with constraints on the den sity of the Local Group 
medium by various stud i es (e.g., iRasmussen fc Pedersenl 
l200ll ; ISembach et~al"1l2003l ; IWilliams et al.ll2005h . ~ 

It is interesting to note that the outer edge of the H I 
disc of NGC 300 appears to be close to balance between 
gravitational and ram pressure forces over a range of rea- 
sonable values for density and velocity. This could indicate 
that the edge of the disc is defined by the aforementioned 
truncation radius at which ram pressure forces start to dom- 
inate over gravity. A stable disc configuration beyond that 
radius would be impossible, suggesting that under certain 
conditions ram pressure could in general be responsible for 



establishing the outer edge of H I discs of galaxies in group 
environments. 

Furthermore, ram pressure could also result in a sig- 
nificant vertical displacement of disc gas and thereby con- 
tribute to the warping observed in the outer disc of many 
spiral galaxies. This situation is shown schematically in the 
left-hand panel of Fig. [15] where the gravitational acceler- 
ation, a z , of the NFW halo perpendicular to the gas disc 
is shown as a function of height, z, above the disc plane. 
Under the assumption of a face-on wind caused by the sur- 
rounding IGM, there will be for each galactocentric radius 
an equilibrium height, z oq , at which the gravitational pres- 
sure in the disc will be in balance with the ram pressure of 
the wind. Of course, there will be no such equilibrium for 
radii greater than the truncation radius beyond which the 
ram pressure always exceeds the gravitational pressure at 
any height above the disc. 

The right-hand panel of Fig. [15] shows the expected 
value of z eq as a function of galactocentric radius in NGC 300 
for two different IGM densities of tt-igm = 10 _J and 
10 -4 cm -3 and under the assumption of a constant rela- 
tive velocity of v Te \ = 150 kms -1 . Over much of the disc of 
NGC 300 the expected equilibrium height, and thus vertical 
displacement, of the gas is very small. Towards the outer 
disc the equilibrium height exceeds 1 kpc with a maximum 
value of z cq w 1.6 kpc in the outermost radial bin under 
the assumption of wigm = 10 -5 cm -3 . If we increase the 
IGM density by a factor of ten to 10 -4 cm -3 the equilib- 
rium height increases to almost 10 kpc at r = 16.6 kpc, and 
beyond that radius there is no solution at all because we 
would exceed the truncation radius. These estimates illus- 
trate that ram pressure may have an impact on the vertical 
structure of the H I of NGC 300 depending on the properties 
of the surrounding IGM. We can also rule out higher IGM 
densities of well above 10 -5 cm -3 as they would result in a 
destruction of the outer parts of the observed H I disc. 

Please note that in all our calculations we neglected 
the gravitational potential of the stellar and gaseous discs, 
hence underestimating the restoring force. This will particu- 
larly be the case for the inner disc of NGC 300 and locations 
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Figure 15. Left: Schematic overview of the gravitational acceleration, a z , of the NFW halo perpendicular to the disc as a function of 
height, z, above the plane. At z = z cq , gravitational pressure and ram pressure are in balance. Right: Height, z eq , above the disc plane of 
NGC 300 where equilibrium between gravitational pressure and ram pressure is achieved. The situation for two different IGM densities 
°f "igm = 10 — 5 and 10 -4 cm -3 is shown under the assumption of face-on motion and a relative velocity of NGC 300 with respect to 
the surrounding IGM of ti re i = 150 kms -1 . 



close to the disc plane, where the ratio between gravitational 
pressure and ram pressure will be significantly higher than 
indicated in Fig. [14] Near the outer edge of the Hi disc, 
however, the contribution from stars and gas to the gravita- 
tional potential is much smaller compared to the dark mat- 
ter component. While the stellar mass surface density across 
the outer disc is immeasurably small, the contribution of the 
gas disc to the gravitational potential is comparable to that 
of the dark matter halo within a few hundred pc of the disc 
plane and negligible for z > 1 kpc. Therefore, our general 
conclusions about the stability of the outer disc are not af- 
fected by the additional contribution from the gas, and our 
numerical results near the outer edge of the disc are accurate 
within a factor of about 2 to 3. 

Our results suggest that the edge of the gas disc of 
NGC 300 may indeed be subject to ram pressure stripping, 
although the effect will likely be subtle. This result is con- 
sistent with our Hi observations which show notable asym- 
metries between the north-western and south-eastern part 
of NGC 300 and signs of distortion along the north-western 
edge of the Hi disc. In this scenario the proper motion of 
NGC 300 on the sky would be in the south-eastern direction, 
pointing away from NGC 55. 

The ram-pressure scenario is also consistent with the 
offsets between the H I gas and the stellar distribution 
found in severa l dwar f galaxies of the Sculptor group by 
I Bouchard et alj (|2005t )Fl Although these effects are by far 
not as strong as in massive galaxy clusters, our results sug- 
gest that even in modest-sized galaxy groups ram pressure 
interaction does play a role in shaping the faint, outer gas 
discs of spiral galaxies in a density regime often not ex- 
posed by shallow Hi observations. This finding highlights 
once more the need for deep and large-scale imaging of galax- 
ies in the 21-cm line of Hi in order to understand the role 
of interaction and feedback between galaxies and the IGM. 



6 It should be noted, though, that their ATCA detections are of 
low significance, resulting in large uncertainties in the position of 
the H I emission. 



5.3 Gas accretion 

Alternatively, some of the observed asymmetries in the H I 
disc of NGC 300 could be the result of cold-mode accre- 
tion flows of primordial gas. Lopsidedness of the stellar 
and gaseous discs o f galaxies is a well-known and common 
phenomenon (e.g.. iBaldwin. Lvnden-Bell fc Sancisll Il98d; 



Richter fc Sancisj|l994l ; lRix fc Zaritskvlll995l; lHavnes et al.l 
199ct ). Among other processes ( Jog fc Combesll2009h . cold- 



mode accretion flows of primordial gas have been dis- 
cussed as the potentia l origin of the obse rved asymmetries 
jBournaud et al.ll2005l ; ISancisi et al.ll2008l). 

Hydrodynamical simulations by iKeres et al l (|2005h 
have shown that there are two distinct modes of gas ac- 
cretion on to galaxies: hot-mode accretion of shock-heated 
gas at temperatures in excess of the virial temperature of 
the dark matter halo, T > T vir w 10 5 to 10 7 K, and cold- 
mode accretion of gas at lower temperatures of T < T v i r . 
While in their simulations the cold mode dominates accre- 
tion at higher redshifts, the hot mode becomes increasingly 
more important at lower redshifts, as many dark matter 
haloes exceed a critical virial mass of about 1Q 11 - 3 - 12 M 
above which accreted gas becomes predominantly shock- 
heated (|Keres et all 120051 . 120081 : iDekel fc Birnboiml 120061 ). 
Furthermore, cold-mode accretion is geometrically different 
from hot-mode accretion in that it occurs along filaments 
and therefore is not radially symmetric. 

Asymmetric cold-mode accretion could have resulted in 
asymmetries in the outer gas disc of NGC 300, because the 
accreted gas would end up in non-circular orbits with sub- 
stantial orbital time-scales in excess of about 1 Ga near the 
outer edge of the gas disc. This scenari o is su pported by nu- 
merical simulations of lBournaud et all (|2005h who conclude 
that accretion along filaments of the cosmic web with accre- 
tion rates of a few Mq a -1 is likely to be responsible for the 
ubiquity of lopsidedness observed in the stellar and gaseous 
discs of galaxies. 

We conclude that sensitive Hi observations of a large 
sample of galaxies would provide the opportunity to sys- 
tematically study the kinematics of gas in the outer disc to 
determine the contribution of external processes such as ram 
pressure or gas accretion. 
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6 SUMMARY 

We have used the Australia Telescope Compact Array to 
map a large area of about 2° x 2° (corresponding to a 
projected size of about 65 x 65 kpc 2 ) around the Sculptor 
group galaxy NGC 300 in the 21-cm line of Hi. We achieved 
a 5cr Hi column density sensitivity of 1.0 x 10 19 cm -2 
over 8 kms - for emission filling the synthesised beam of 
180" x 88" FWHM. The corresponding 5a Hi mass sensi- 
tivity is 1.2 x 10 5 M under the assumption of a distance of 
1.9 Mpc for NGC 300. 

The major results and conclusions from our analysis of 
the Hi emission in NGC 300 are: 

• NGC 300 is characterised by a dense inner disc with 
column densities in excess of about 5 x 10 20 cm -2 and an 
extended outer disc with column densities of typically less 
than 2 x 10 20 cm -2 . The inner disc has an orientation angle 
of 290° and is perfectly aligned with the optical (stellar) 
disc. The outer disc has a totally different orientation angle 
of 332°, and its major axis spans more than 1° on the sky 
equivalent to a linear diameter of about 35 kpc within the 
10 19 cm" 2 column density contour. 

• Under the assumption that the gas is optically thin we 
derive a total H I mass for NGC 300 of 1.5 x 10 9 M Q . This is 
only slig htly short of the single- dish value of 1.7 x 10 9 Mq 
found bv lKoribalski et al.1 (|2004h . indicating that the ATCA 
recovered about 90 per cent of the total flux. Of the total 
Hi mass of NGC 300, about 8.5 x 10 s M Q (or 57 per cent) 
is contributed by the inner disc and 6.5 x 10 8 Mq (or 43 per 
cent) by the outer disc. 

• We fitted a tilted ring model to the velocity field of 
NGC 300 to determine the Hi rotation curve out to a ra- 
dius of 18.4 kpc. The derived rotation curve rises out to a 
radius of approximately 10 kpc where we observe a maxi- 
mum velocity of about 100 kms -1 . Further out the rotation 
curve slowly decreases to about 83 kms -1 in the outermost 
ring at 18.4 kpc radius. 

• We constructed different mass models of NGC 300 in- 
volving a stellar disc, a gas disc, and a dark matter halo. 
We obtain good fits for both the Burkert and NFW dark 
matter halo models, both of which successfully cope with 
the decreasing rotation velocity in the outer parts of the 
galaxy. Fits with either Newtonian dynamics or Modified 
Newtonian Dynamics (MOND) and no dark matter halo at 
all result in much higher Xicd values and unrealistic mass 
scaling factors for the gaseous disc. 

• The twisted structure of the H I disc of NGC 300 sug- 
gests distortion by tidal interaction, but no nearby compan- 
ion galaxy is currently known that could have exerted strong 
tid al forces on the disc. At the same time, optical studies 
by iBland-Hawthorn et al.l l|2005h provide evidence against 
tidal interaction by demonstrating that the stellar disc of 
NGC 300 has a simple exponential radial profile out to the 
largest detected radii without any sign of distortion. 

• Optical observat io ns by 

IVlaiic. Bland-Hawthorn fc Freeman] (|2009h have shown 
that the stellar population in the outer disc is predom- 
inantly old (> 1 Ga) and that therefore the outer gas 
disc does not currently play a significant role in the star 
formation activity of NGC 300. 

• Asymmetries in the outer H I disc suggest that NGC 300 
is affected by ram pressure interaction while moving through 



the intergalactic medium of the Sculptor group, with its 
proper motion on the sky directed approximately to the 
south-east. Estimates show that, under reasonable assump- 
tions on the IGM density and relative velocity, the outer 
edge of the Hi disc of NGC 300 may be significantly dis- 
torted by ram pressure. From our analysis we conclude that 
the density of the IGM in the Sculptor group must be lower 
than about 10 -5 cm -3 as otherwise the outer parts of the 
H I disc of NGC 300 could not have survived. 
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